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Using (106 ±4) x 10 6 0(3686) events accumulated with the BESIII detector at the BEPCII e + e" 
collider, we present the first measurement of decays of \a to vector meson pairs 00, ojuj, and uxj>. 
The branching fractions are measured to be (4.4 ± 0.3 ± 0.5) x 10~ 4 , (6.0 ± 0.3 ± 0.7) x 10~ 4 , 
and (2.2 ± 0.6 ± 0.2) x 10~ 5 , for Xci <t>4>, ojuj, and uj(j), respectively, which indicates that the 
hadron helicity selection rule is significantly violated in Xcj decays. In addition, the measurement 
of XcJ — ► w0 provides the first indication of the rate of doubly OZI-suppressed \ c j decay. Finally, 
we present improved measurements for the branching fractions of XcO and Xc2 to vector meson pairs. 

PACS numbers: 14.40. Pq, 12.38.Qk, 13.20.Gd, 13.25.Gv 
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Decays of the \c,j (J — 0,1,2) P-wave charmonium 
states are considered to be an ideal laboratory to test 
QCD theory. The initial theoretical calculations of \ c j 
exclusive decays into light hadrons predicted branching 
fractions that were smaller than the experimental mea- 
surements [l[ . With the inclusion of the color-octet mech- 
anism 0, calculations of \cj decays into pairs of pseu- 
doscalar mesons and pairs of baryons came into reason- 
able agreement with the experimental measurements, in- 
dicating the importance of the color-octet mechanism. 

In the case of XcJ decays into pairs of vector (J PC = 
1 ) mesons VV, where V is an uj or 0, the branch- 
ing fractions for XcO/2 decays to 00 and have been 
measured to be at the 10~ 3 level [!, 0|, which is much 
larger than predictions based on perturbative QCD calcu- 
lations [H . Decays of the Xci into 00, ujuj and ujcj) violate 
the helicity selection rule (HSR) and are expected to be 
highly suppressed @. In addition, the decays XcJ — > ^<P 
are doubly OZI suppressed and have yet to be observed. 
Recently, long-distance effects in Xci decays 0, H[ have 
been proposed to account for the HSR violation. Precise 
measurements of Xci ~ > W decays will help clarify the 
influence of long distance effects in this energy region. 

In this Letter, we report measurements of XcJ decays 
into 00, ujuj, and uj<f> modes, where is reconstructed 
from K + K~ or 7r + 7r~7T°, u> from 7T + 7r~7r°, and 7T° from 
77. The data samples used in this analysis consist of 
(106±4) x 10 6 ip' decays and 42.6 pb _1 of continuum data 
at y/s = 3.65 GeV acquired with the BESIII detector 



The cylindrical core of the BESIII detector consists of a 
helium-gas-based Main Drift Chamber (MDC), a plastic 
scintillator Time-of-Flight system (TOF), a CsI(Tl) Elec- 
tromagnetic Calorimeter (EMC), and a muon counter. 
The charged particle and photon acceptance is 93% of 
47r, and the charged particle momentum and photon en- 
ergy resolutions at 1 GeV are 0.5% and 2.5%, respec- 
tively. The BESIII detector is modeled with a Monte 
Carlo (MC) simulation based on GEANT4[lol. The 
optimization of the event selection and the estimation 
of physics backgrounds are performed with Monte Carlo 
simulations of '0(3686) inclusive/exclusive decays 

The final states of interest are 72(i^ + i^~), 572(7r + 7r - ), 
and 3"/K + K ~tt + tt~. Event candidates are required to 
have four well reconstructed charged tracks with net 
charge zero, and at least one, five, or three good pho- 
tons, for 00, ujuj , and w0, respectively. 

Electromagnetic showers in BESIII detector are recon- 
structed from clusters of energy deposits in the EMC. 
The energy deposited in nearby TOF counters is included 
to improve the reconstruction efficiency. A good pho- 
ton is a shower in the barrel region (| cos6*| < 0.8) with 
at least 25 MeV energy deposition, or in the endcaps 
(0.86 < I cos6»| < 0.92) with at least 50 MeV energy de- 
position, where 6 is the polar angle of the shower. Show- 
ers in the region between the barrel and the endcaps are 
poorly measured and excluded. Timing requirements are 
used in the EMC to suppress electronic noise and energy 
deposits unrelated to the event. 
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Charged tracks are reconstructed from MDC hits. 
Each charged track is required to be in the polar an- 
gle region |cos#| < 0.93 and to pass within ±10 cm of 
the interaction point in the beam direction and within 
±1 cm in the plane perpendicular to the beam. 

A kinematic fit constrained by the initial e + e~ four- 
momentum in the lab frame is applied to the decay 
hypotheses -0(3686) -> j2(K + K~), 5j2(tt + tt~), and 
3 r yK + K^TT + TT^ . The final state photons are identified 
with the photon-charged-track combination that has a 
minimum x\c va hie (for definition of x\c j see ) wnen 
sampling all candidate photons. The vertex of all charged 
tracks must be consistent with the measured beam in- 
teraction point. The x\c selection efficiency is opti- 
mized using the ratio of signal to backgrounds in the 
data: xlc < 60 for 1^{K + K~), 3jK+K-ir+ir- , and 
x\q < 200 for 572(7r + 7r - ) is required. To separate the 
K* 1 from 7r~^~ in the 3~/K + K~7r + Tr~ final state, two kaons 
are identified with the requirements that P(K) > P(n) 
and P(K) > P(p), where P(X) is the probability of hy- 
pothesis X as evaluated from the TOF and dE/dx infor- 
mation. 

The mass windows for resonance candidates are set ac- 
cording to the optimized ratio of signals to backgrounds 
in the data. The tt° candidates are selected by requir- 
ing 0.1 < M 77 < 0.15 GcV/c 2 . The <j> and uj can- 



didates are selected by requiring \Mk+k- — 1-019| < 
0.015 GeV/c 2 , \M„+„-^o - 1.019| < 0.030 GeV/c 2 , and 
\M^ +w - %0 - 0.783| < 0.050 GeV/c 2 , for -> K+K~, 
(f) — > 7T + 7r _ 7r°, and uj — > ir + ir~ 

For XcJ-KM>->2(K + K~'. 
the minimum value of (M (1) 



k , respectively, 
the two </> candidates with 

K+K _-i.mf + {Mp +K _- 

1.019) 2 are taken as the signal. No artificial 0-pair peaks 
are produced when this selection criteria is applied to 
MC simulation of the process Xcj —> 2(K + K~). A scat- 
tcrplot of masses for one K + K~ pair versus the other 
K + K~ pair is shown in Fig. [TJa), where a clear <jxp sig- 
nal can be seen. The Mk+k- distribution, after requir- 
ing that the other two kaons are consistent with being 
a </>, is shown in Fig. [TJb). A <fi peak is clearly seen 
with very low background. The <p<p invariant mass dis- 
tribution for the selected events is shown in Fig. [5] (a), 
where XcJ signals are clearly observed. The MC simu- 
lation shows that the peaking backgrounds, i.e., back- 
grounds that produce Xcj signal peaks, are mostly from 
XcJ 4>K + K~ and 2{K + K~) final states; the back- 
grounds from misidentified charged particles are negligi- 
ble. The levels of the peaking backgrounds are evaluated 
from N AB = r A Nf-r B N%, where Nf{N B t ) is the num- 
ber of data events falling into box A (B), as indicated in 
Fig. [TJa), and the normalizing factors r, = N^f g /N^ c 
with i = A or B are determined from MC simulation 
for modes XcJ — > 4>K + K~ and 2(K + K~), respectively. 
Here N^ g c (N^°) is the number of MC events falling 
into the signal box (A or B). These backgrounds will 
be indistinguishable from signal events; therefore, we fix 
their normalization, independently for each Xcj peak, in 
the final fit. 




(e) M o (GeV / c ) 



(f) M „ (GeV / c ) 



FIG. 1: The left column shows scatterplots for events within 
the XcJ mass region. The boxes indicate the signal region 
(without label) and sideband regions labeled as A and B. The 
plots in the right column are the one-dimensional projections 
of the system recoiling against a selected <j> or u> resonance. 
Plots (a) and (b) are for the y2(K + K~) mode; (c) and (d) for 
the 572(7r + 7T~) mode; and (e) and (f) for the 3~/K + K~tt + tv~ 
mode. 



To study Xcj — > ww decays into the 2(7r+7r n°) final 
state, two 7T° candidates are selected by minimizing the 
value of (M$ - 0.135) 2 + (M$ - 0.135) 2 when sam- 
pling all four-photon combinations from the selected five 
photons. The 7r + 7r _ 7r° combination closest to the nom- 
inal uj mass is taken as one uj candidate, and the re- 
maining three pions are assumed to be from the other 
uj. No artificial cj-pair peaks are produced from the 
application of this w-selection criteria to a MC simula- 
tion for Xcj — > 2(7T + 7r~7r°). A scatterplot of the mass 
for one 7r + 7r~7r° pair versus the other 7r + 7r~7r° pair is 
shown in Fig. [He), and the M^+^-^o distribution for the 
three pions recoiling against an u> candidate is plotted 
in Fig. [Ud). The ujuj mass spectrum is shown in Fig. [2] 
(c), where XcJ signals are prominent. The MC simu- 
lation shows that the backgrounds in the ujuj signal re- 
gion include peaking backgrounds from Xcj — > Luir + ir~7r a 
and 2(7r + 7T _ 7r°), and non-peaking backgrounds from the 
0(3686) decays into the same final states without inter- 
mediate Xcj states. The backgrounds from misidentified 
charged particles are negligible. Potential backgrounds 
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FIG. 2: Invariant mass of VV for (a) (j>4> mode in the 
y2(K + K~) final state, (b) 4>4> mode in the -y-ir + Tv~ tt° K + K~ 
final state, (c) ujuj mode in the r y2(ir + tt~ ir°) final state, and 
(d) co(f> mode in the ^tt + tt~tt°K + K~ final state. The points 
with error bars are the data; the solid lines are the fit re- 
sults; and dotted lines represent the signal components. The 
shaded and open histograms in (a,b) and (c) , respectively, are 
peaking backgrounds. In (c), the shaded histogram denotes 
the non-Xcj backgrounds. In (d) the long dash line is back- 
ground normalized by a simultaneous fit to iocp sidebands, and 
the dash-dot line is non-^cj background. 



from Xcj -> # -> 2(7r+7T tt°) and x c o/2 -> VV ~> 
2(7r + 7T _ 7r°) do not survive our selection criteria. As in 
the Xcj mode, the sizes of the peaking backgrounds 

from Xcj — > Luir + ir~ir a and 2(7r + 7r _ 7r°) are evaluated by 
selecting data events located in sideband boxes A and 
B, respectively, as indicated in Fig. [ljc). The peaking 
backgrounds are normalized according to the ratio of MC 
events falling into the signal region and those falling into 
the sidebands. The normalization of these peaking back- 
grounds is fixed in the final fit. 

To study Xcj u<fi and <fxj> decays into the 

K + K~ir + ir~ir° final state, the photon pair with invari- 
ant mass closest to the ir° nominal mass is taken as the tt° 
candidate. A scatterplot of masses for K + K~ pairs ver- 
sus that for 7r + 7r~7r° pairs is shown in Fig.QJc), and the 



M- 



distribution for events satisfying 



K + K~ 



is shown in Fig. [TJf), where the ui — > 7r + 7r~7r° and 
cf) — > 7r + 7r~7T° signals are clearly seen. The (jxfi and ui<j> 
mass spectra are shown in Figs. [2] (b) and|2](d), respec- 
tively. Similar to the case for Xcj 2(K + K~), 
the peaking backgrounds from the Xcj 07r + 7r~7r° or 
4>K + K~ , and K + K~tt + tt~tt° are evaluated by selecting 
data events falling into sideband boxes A and B, respec- 



tively, as indicated in the inserted plot in Fig. QJe) . The 
peaking backgrounds are normalized according to the ra- 
tio of MC events falling into the signal region and those 
falling into the sidebands. The normalization of these 
peaking backgrounds is fixed in the final fit. 

The numbers of observed events are obtained by fit- 
ting the Myv distributions. The observed line shapes 
are described with modified Xcj MC shapes plus back- 
grounds. Possible interference effects between the sig- 
nal mode and the peaking background modes are not 
considered for all modes. The original Xcj MC shapes 
are generated by a relativistic Brcit-Wigncr incorporated 
with full hclicity amplitudes in the EvtGen package 14| , 
and their masses and widths are set to the nominal val- 
ues [13 ■ In the fits they are modified by convolving 
them with Gaussian functions G(Myv — SMj, crj), where 
5Mj and aj correct the Xcj mass and width or res- 
olution, respectively, in the simulation. The values of 
5Mj and crj, determined from the fits, are less than 1 
McV for all modes and from 1 to 5 MeV, respectively. 
Backgrounds from QED processes, which are estimated 
from the application of a similar analysis to the contin- 
uum data, are negligible. For Xcj — > 4>4>i t ne peaking 
backgrounds are fixed to the sideband estimates as men- 
tioned above, and other combinatorial backgrounds are 
parameterized by a second-order polynomial with param- 
eters that are allowed to float in the fit. For all modes, 
a maximum-likelihood technique [lii] is employed to es- 
timate parameters. After projecting the best fit into 
the binned histograms shown in Fig. ® we determine 
x 2 /ndf = 0.46 for X cj -> H -> 2{K + K~) and 0.50 for 
the XcJ ~^ 4>4> ~~ ^ K + K~7t + tt~tt , where ndf is the num- 
ber of degrees of freedom. The fitted results are plotted 
in Fig.[3Ja) and (b), respectively. The numbers of signal 
events are listed in Table fl] 

For the XcJ ~ * channel, backgrounds include the 
peaking backgrounds estimated from uj sidebands indi- 
cated in Fig. [T] (c), non-^cj backgrounds (-0(3686) — > 
^louj) fixed at the normalized MC shape of phase space 
using the data information, and smooth combinatorial 
backgrounds that are parametrized by a second-order 
polynomial. The x 2 / n df for the fit is 0.97. The fit results 
are shown in Fig. [2] (c). 

To extract the signal yield, as well as to estimate the 
statistical significance for the Xcj —> mode, a si- 
multaneous fit is performed to M U! f, distributions both 
in uxp signal and sideband regions of boxes A and B 
[see Fig. Q] (c)]. The peaking backgrounds are normal- 
ized according to the ratio of MC events falling into 
the signal region to those falling into the sideband re- 
gions for the -0(3686) —> j(f>ir + ir~ir a , jujK + K~ and 
■0(3686) — > ~fK + K~ir + Tr~ir° events that are within the 
Xcj mass region. Because of the low signal yield in this 
mode, the parameters 8Mj and a,j of the modified MC 
shapes are fixed at the values determined in the fit of 
Xc.j -> H -> K+K-tt+tt-tt . The x 2 /ndf is 0.62. The 
fit results are shown in Fig. [2] (d), and the numbers of 
signal events are listed in Table |U 
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TABLE I: Summary of the branching fractions (B) for Xcj — > 
(jxj), ujui, and oj<f>. Here -/V ne t is the number of signal events, 
e is the detection efficiency. The upper limit is estimated at 
the 90% C.L. 
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9.5 + 0.3 + 1.1 


Xd — > 


597 ± 29 


13.2 


6.0 + 0.3 + 0.7 


Xc2 ^ UJUJ 


762 ± 31 


11.9 


8.9 + 0.3 + 1.1 


-> 2(7r+7r-7r°) 








XcO -* W0 


76 ±11 


14.7 


1.2 ± 0.1 ± 0.2 


Xd — > U(f> 


15 ±4 


16.2 


0.22 + 0.06 + 0.02 


Xc2 — > LO(j) 


< 13 


15.7 


< 0.2 











The uncertainties due to the modified \cJ MC shapes 
are estimated by replacing them with Brcit-Wigncr func- 
tions convolved with the instrumental resolution func- 
tions in the fits. The quality of the resulting fit is not as 
good as using the modified MC shapes. The difference of 
signal yields varies from 1% to 4%, and this is included 
as a systematic error. 

The detection efficiencies are determined from MC sim- 
ulations for the sequential decays -0(3686) — > 7Xcj — > 
VV, V decays into the selected final state. The decays 
■0(3686) — > 7XcJ are generated by assuming a pure El 
transition. The \ c j — > VV decays and subsequent de- 
cays of the V are modeled with helicity amplitudes that 
provide angular distributions consistent with the data. 

The systematic uncertainties on the Xcj decay branch- 
ing fractions arise from the ir^ and tracking, K 
identification, EMC shower reconstruction, number of 
-0(3686) decays, kinematic fitting, modified MC shapes, 
background estimation, Xcj signal extraction and uncer- 
tainties from branching fractions of -0(3686) — > ^Xcj, 
— > K + K~, lu — > 7r + 7r~7T° and tt° — > 77. The un- 
certainties caused by MDC tracking are estimated to be 
2% for each charged track [l7|. The uncertainty due 
to identification is evaluated to be 2% per kaon 
The uncertainty due to the photon reconstruction 
is determined to be 1% for each photon (T3). The un- 
certainty in the number of 0(3686) decays is 4% [i"2| . 
The uncertainties due to the kinematic fit are deter- 
mined by comparing the efficiency at the given x\c va ^" 
ues for the MC sample to control samples selected from 



data, i.e, ^(3686) -> 7 00 -> j2(K + K~), 0(3686) -> 
ir TT°J/iP, J/4' -> 2(7r+7r-),7r°2(7r+7r-) and 0(3686) -> 
7T + 7r _ J/tp, J/i/j —> K + K~tv°. The kinematic-fit un- 
certainty varies from 0.5% (72(7r + 7r _ 7r°) mode) to 3.7% 
("fK + K~ir + 7T~TT Q mode). The uncertainties of the peak- 
ing backgrounds for Xcj ~> 00 — * 2{K + K~) are evalu- 
ated by comparing the sideband estimates to the exclu- 
sive MC simulation on the modes Xcj — > <t>K + K~ and 
2(K + K~). while for other modes the uncertainties are 
estimated by varying the size of sideband boxes. The un- 
certainties of the peaking background estimates are less 
than 3%. The uncertainty from the MC normalization 
factor is found to be negligibly small. The total system- 
atic uncertainties arc 10% for X cj — > 00 — > 2(K + K~) 
mode, and 11% for Xcj —> uuj —> 2(-7r + 7r~7r ), Xcj — ► 
00, Ld(j) —> K + K~ir + 7T~ir modes. 

The branching fractions for Xcj decays are determined 
from B = iV nc t/(A^(3686)erL^); where N nQt and e are 
the number of net signal events and the detection effi- 
ciency, respectively. The detection efficiencies are listed 
in Table ffl Here ^ (3686) = (106 ± 4) x 10 6 \V2} is 
the number of 0)(3686) events, and Yii &i is the prod- 
uct of world average branching fractions values [15J for 
■0'(3686) — > jXcj and the other meson decays that are 
involved. For the Xcj — > 00 — > K + K~7t + tt~it branch- 
ing fraction we double the efficiency listed in Table|T]since 
our analysis sums over the two combinations for each to 
decay to either K + K~ or ir + ir~ir . The resulting branch- 
ing fractions are listed in Tabic HI The statistical signifi- 
cance of Xcj — > ^0 is derived from the change of — 2 ln£ 
obtained from fits with and without each of the three 
Xc.r ^0 signal components. Wc obtain a significance 
of 4.1cr for Xci - ► ^0 and 1. 5a for Xc2 u>4>- The signif- 
icance of the XcO ~~ ^ w signal is 10er. Using the Bayesian 
method, the upper limit for the number of signal events 
of the Xc2 +> w0 mode is 13 at the 90% confidence level 
(C.L.). The branching fractions for \ c j cjxj) mea- 
sured in 2{K+K~) and (A'+A'-)(7r+7r~7r°) final states 
are combined into a weighted average, where common 
systematic uncertainties are counted only once. 

In summary, the HSR suppressed decays of Xci — > 
00, ujlo. and the doubly OZI suppressed decay XcO — > ^4> 
are observed for the first time. The branching fractions 
arc measured to be (4.4 ± 0.3 ± 0.5) x 10" 4 , (6.0 ± 0.3 ± 
0.7) x 10~ 4 , and (1.2 ± 0.1 ± 0.2) x 10" 4 , for X a -> <t><P, 
loco 7 and XcO — > UJ( f>i respectively, We also find evidence 
for Xci — > ^0 decay with a signal significance of 4.1er. 
The branching fractions for XcO/2 ~^ 00, decays are 
remeasured with a precision that is better than those 
of the current world average values [l5| . These precise 
measurements will be helpful for understanding Xcj de- 
cay mechanisms. In particular, the measured branching 
fractions for Xci — > VV indicate that HSR is significantly 
violated and that long distance effects play an important 
role in this energy region. The long distance effects from 
the intermediate charmed meson loops in Xci — > 4>4> an d 
uju) decays @,Q can contribute to the branching fractions 
at the level of 10~ 4 but are more than an order of mag- 
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nitude too small to explain the doubly OZI suppressed 
decay rate for Xd ~^ that we measure ||. 
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